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Abstract 



We investigate linear parabolic maps on the torus. In a generic case these maps 
are non-invertible and discontinuous. Although the metric entropy of these systems 
is equal to zero, their dynamics is non-trivial due to folding of the image of the unit 
square into the torus. We study the structure of the maximal invariant set, and in a 
generic case we prove the sensitive dependence on the initial conditions. We study the 
decay of correlations and the diffusion in the corresponding system on the plane. We 
also demonstrate how the rationality of the real numbers defining the map influences 
the dynamical properties of the system. 
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1 Introduction 

Linear area-preserving maps on the torus T 2 are often analyzed in the theory of ergodic 
Hamiltonian systems jl], 01. The map M, defined through 



ax + by | m od 1 
cx + dy \ mod 1, 



can be represented by the Jacobian matrix 



a b 
c d 



(1) 



(2) 



Due to the area-preserving condition ad — be = 1 the matrix g pertains to the non- 
commutative group SX(2, R,). Let us emphasize here that the area-preserving property 
holds for the map when unfolded on the plane. If the dynamics takes place on the torus, 
folding causes overlap of some parts of the image of the unit square, so the area is not 
conserved. 

For each matrix g there exist its inverse g^ 1 = [d, —b; — c, a], which defines the inverse 
map M when the dynamics takes place on the plane (x,y). However, if the dynamics 
takes place on the torus T 2 (modulo 1 restriction in Eq. (E)), one has to wrap the image 
of the basic square back into the torus, and some parts of this image may overlap. Conse- 
quently, some points in the torus may not have any pre-images, and the map ([!]) may not 
be invertible on the torus. On the other hand, the systems defined by matrices with all 
integer elements are invertible. The set of these matrices Gi n t forms a discrete subgroup of 
SL(2, R). Also matrices of the type [1, t\; 0, 1] and [1, 0; t2, 0] represent invertible (but not 
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continuous) maps on the torus for any non-integer parameters t\ and ti- These matrices 
correspond to the horocyclic flows (lj and form two continuous subgroups of SL(2, R), 
denoted by Gjj x and Gh 2 - 

Dynamical properties of the map (|l|) can be characterized by the trace T =Tv(g) = 
a + d. For \T\ > 2 the map is hyperbolic, eigenvalues of g are real; z > 1 and 1/z (see e.g. 
H ) . There exist infinitely many periodic orbits and all of them are unstable with the same 
stability exponent Ai = ln(z). The Arnold cat map defined by the matrix g = [1, 1; 1, 2] is 
one of the most celebrated examples of such chaotic dynamical system. Since all elements 
of g are integer, this map is invertible on the torus T 2 (it is an automorphism on the 
torus). For \T\ < 2 the map is called elliptic. The eigenvalues of g are complex and 
both Lyapunov exponents are equal to zero. Dynamics of such a system is regular and 
corresponds to a rotation. Properties of the elliptic maps on the torus and some features 
of the elliptic-hyperbolic transition were studied by Amadasi and Casartelli H 

The intermediate case, \T\ = 2 is called parabolic. Although this case is called 
by Berry || 'special, non-generic, set-of-measure-zero, infinitely-improbable-unless-you- 
deliberately-set-out-to-create-them' case, we believe that it is worth analyzing, for two 
reasons. On one hand, these maps are interesting, as they lie in between the elliptic and 
the hyperbolic cases, which are very important for several physical applications. On the 
other, the parabolic maps display several remarkable properties. In particular, we show 
that a generic linear, area-preserving, parabolic map on the torus displays sensitivity on 
initial conditions. Furthermore, we demonstrate how the rationality of numbers defining 
the map affects the dynamics of the system. A related study of the systems leading to the 
interval exchange maps were discussed in j|. 

This paper is organized as follows. In section 2, we discuss the general properties of 
parabolic maps on the torus and show corresponding families of ID maps. In section 
3 and 4, we analyze the case of rational and irrational maps, respectively. In section 
5, we analyze the property of sensitive dependence on the initial conditions. Decay of 
correlations and the diffusion rate are analyzed in section 6. 

2 Parabolic maps on the torus 

If the trace T of the matrix (H) fulfills |T| = 2, the linear map on the torus T 2 is called 
paraboIic.[| The corresponding dynamics is not chaotic and describes a shear flow. Since 
trace of a product of two matrices is usually not equal to the product of their traces, the 
set C?2 of all matrices g with the trace equal to 2 and the determinant equal to unity (or 
equivalently, the set of all area-preserving, linear parabolic maps on the torus) does not 
possess a group structure. Apart from the matrices with integer entries and trace equal 
to two (such as [3, — 2; 2, — 1]), the matrices belonging to the groups Gh-i and Gh 2 a ls° 
correspond to invertible parabolic maps. If both elements b and c of the matrix (^) are 

1 To avoid the confusion let us mention that the 2D maps analyzed in this work are entirely different 
from ID maps constructed out of parabola, also referred in the literature as parabolic maps. 
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non-zero, any matrix g £ G2 may be represented as 



9 



l + A A/a 
-Aa I- A 



(3) 



/ _ „ m _i „.u , A„, (5) 



where both the free parameters A and a are real. Both eigenvalues of g are equal to unity, 
and the eigenvector reads v = (1, —a). For parabolic maps, the two dimensional dynamics 
on the torus can be decomposed into a family of ID maps, which are defined on the lines 
parallel to the vector v. The character of the number a plays therefore a special role: for 
rational a = r/s, any line parallel to v forms a closed circle on the torus, consisting of at 
most r + s lines on the unit square. On the other hand, for any irrational value of a, any 
line parallel to v winds densely the entire torus. 

To reduce the 2D dynamics into a family of ID maps, let us introduce new variables 
{p, q} defined as 

h = x ( 4) 

y q = ax + y, 

where p £ [0, 1), and q E [0, 1 + a) . In these variables, the map M takes the following 
form: 

p' = p + Aq/a I mod 1 

[{1 - A)q - ap] - a[p + £q], 

where \p] denotes the integer part of the real number p. This formula is derived as 
following. Consider the line joining point (x,y) and its image (x',y') in the xy-plane. It 
crosses the horizontal discontinuities (at integer y values) [y'\ times, and the vertical ones 
[x'] times. Then the above formula is obtained from q' = q — [y'\ — a[x'] since the slope 
of the line is —a. Note that if qo = axo + yo is the initial q value, the variable q can be 
represented as q = qo + l\ + /2 a with integers l\ and I2, and the dynamics of q can be 
regarded as a motion on the integer lattice (^i,^)- For some purposes, it is convenient to 
ignore the operation modulo 1 in Eq. (||), which corresponds to considering the dynamics 
on the xy-plane instead of the torus. The dynamics is described by the family of ID maps 

P' = P+^-{qo-a[P}-[q -aP}}, (6) 

where P £ H and the initial condition axo + yo = Qo £ [0, 1 + a) is treated as a parameter. 
The slope of the above map is constant and equal to unity. Thus the Lyapunov exponent is 
equal to zero, in agreement with the assumed parabolicity of the 2D map M. For rational 
a = r/s, the map (g) is periodic with period s. In the simple case a = 1, the map (||) 
becomes a lift of degree one, often analyzed in the theory of dynamical systems (see e.g. 



3 Rational parabolic maps on the torus 

We shall call a parabolic map on the torus rational, if the parameter a, which enters 
the matrix g given by (||), is rational. If a, A and A/ a are rational but not integer, the 
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corresponding map M is discontinuous on all sides of the square and is not invertible on 
the torus. A simple example is given by a = 1 and A = 1/2, which represents the map 
M s ' ^ 

U 2^~^2^l mod 1, 

associated to the matrix g s := [3/2, 1/2; —1/2, 1/2]. Figure ||(a) represents the first iterate 
of the basic square on the plane. Due to the wrapping of this figure back into the torus, 
some parts of it overlap, and some fragments of the torus (denoted in white in Fig. |l|(b) 
do not have any pre-images with respect to M s . Subsequent images of this set are spread 
over entire torus, which is what makes the dynamics of the system interesting. 

Since a = 1, any trajectory originating from a single point is restricted to a circle 
on the torus represented by at most two parallel lines (of the same slope equal to — 1) 
in the unit square. Therefore the map M s does not have an attractor in the sense of 
attracting single trajectories. Figure |] presents a graph obtained from ten thousand 
iterations of trajectories originating from ten thousand different initial points randomly 
chosen according to the uniform distribution on the torus. To avoid transient effects, first 
hundred points of each trajectory are not marked. Denote the maximal invariant set JjJ 
(black in the figures) by 

S= p| M n {T 2 ). (8) 

neN 

This set forms a support of a 'natural' invariant measure of the map M, which may 
be generated by iterating the initially uniform (Lebesgue) measure by the associated 
Frobenius-Perron operator. The measure /i* need not cover the set S uniformly. 

Remark 1. Observe the rotational symmetry of S around the central point Q = 
(1/2, 1/2). Moreover, for any fragment of the set S (black in the figure), one can find the 
corresponding white fragment belonging to S = T 2 \ S. For example, the two large black 
fragments at the top and the bottom of the square form a parallelogram on the torus, 
which corresponds to the white parallelogram consisting of the two large white fragments 
on the left and on the right of the unit square. Thus there exists a symmetry in shifting 
the figure by vector (1/2, 1/2) and interchanging the colors. Therefore, the volume of the 
invariant set S is equal to 1/2. In other words, for a = 1, A = 1/2, the invariant set of 
the ID map (|6]) (see Fig. ||) has volume 1/2 independently of the parameter q. 

The fact that the invariant set S contains a ball of a positive radius implies that its 
dimension is equal to 2. Figure |2](b) presents the magnification of the invariant set in the 
vicinity of the point Q. Note an infinite sequence of interlacing black and white parts 
of the figure which accumulate on the diagonal line and form some kind of self-similar 
structure. 

Remark 2. Numerical experiments suggest that the invariant set S is structurally 
stable. It is thus possible to get a practical approximation of the set S by iterating a single 
trajectory of an modified system close to the analyzed map M s . Perturbing any of four 
elements of g s (say element a by putting a — > a + e with a real parameter < e « 1), 
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we get a perturbed map M £ . Formally, this matrix does not belong to G2, since the area- 
preserving property and parabolicity, T = 2, are e-violated. The perturbation couples all 
ID maps together, so that by iterating one single trajectory we obtain the picture of the 
invariant set S e of the map M e . The numerical results show that this set approximates 
the set S in the Hausdorff metric, while the quality of this approximation improves when 
decreasing the perturbation strength e. 

4 Irrational parabolic maps on the torus 

A generic matrix belonging to G2 contains irrational elements. To analyze such a case in 
some details we take the golden mean 7 := (\/5 — l)/2 for the parameter a entering the 
matrix dH), and set A = 1/2 arbitrarily. This choice corresponds to the map M 7 given by 

X = |x + I mod 1 

V '2% 2^ I mod T 

Due to the irrationality of the parameter a, the line parallel to the eigenvector of g 7 winds 
densely around the entire torus. In other words, the ID map @ is not periodic. Thus, by 
iterating (almost every) single trajectory by the corresponding map M 7 given by Eq. (|l|), 
we obtain full information concerning the invariant attracting set S~. Such a set is shown 
in Fig. 0(a). Numerical analysis allows us to estimate the capacity dimension of this set 
Dq ~ 1.71 ± 0.05 and the correlation dimension D2 ~ 1.70 ± 0.02. These data suggest 
strongly a fractal structure of the set 5 7 . Note that the transformation M 7 is volume- 
preserving on the plane, and the existence of the attractor for this map is solely due to 
the modulo 1 condition, which corresponds to wrapping the image of the unit square back 
on the torus. 

Every irrational number can be approximated by rational numbers by truncating its ex- 
pansion into the continuous fraction. For the golden mean, the consecutive approximations 
are given by the ratios of the adjacent Fibonacci numbers. Taking such approximations 
as 71 = 2/3 and 72 = 5/8, we get two rational maps (g) defined by a\ = 2/3, A\ = 1/2 
and a.2 = 5/8, A2 = 1/2. Estimation of the volume of the corresponding invariant sets S\ 
and S2 gives V\ = 0.3 and V2 = 0.2, respectively. The volume decreases for rational maps 
obtained for better rational approximations of 7 characterized by larger denominators. 
These results suggest that the dimensions of the sets Si and S2 are equal to 2. However, 
an unusually slow convergence of the standard numerical procedure made it difficult to get 
a reliable numerical estimations of these dimensions. Observe similarities between Si and 
S2 (depicted in Fig. ||(c) and (b)) and the invariant set 5 7 of the irrational map presented 
in Fig. §(a). Even though the sets S2 and S~ look similar, we conjecture that their struc- 
ture is very different. Invariant set £2, corresponding to the rational map, is conjectured 
to have a non-zero volume and the capacity dimension equal to 2. This contrasts the 
properties of the set 5 7 invariant under the irrational map M 7 . 
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5 Discontinuity and sensitivity on initial conditions 



Let us denote by D the boundary of the unit square. The linear, parabolic linear area- 
preserving map (|l]) is discontinuous at D with probability one with respect to the Haar 
measure on SL(2, R) (the only case that it is not discontinuous is when all the entries of 
the matrix are integer). The discontinuity can be measured by A, the minimal size of the 
distances d(M(0, y), M(l, y)) and d(M(x, 0), M(x, 1)) on the torus T 2 , which depends only 
on the parameters of a parabolic map. It follows from the definition of A that the condition 
A = is satisfied only by countable number of pairs (A, a). Thus, a generic parabolic 
map is characterized by a positive A. The lines of discontinuity have implications for 
the dynamics of the system, since they split a bundle of neighboring trajectories. Similar 
effects were reported for elliptic maps on the torus in §| @- 



Let us consider the set of all pre-images of the set D. Although the inverse map M 1 
may not exist, we define the pre-images M~ k (D) by 

M~ k (D) = {(x, y) : M k (x, y) G D but M n (x, y) £ D for < n < k}. (10) 

In an analogous way, M k (D) denotes the /c-th image of D with respect to M. Moreover, 
following ||, we define 

D~ = (J M- k (D) and D ± = (J M k (D). (11) 

keN k&Z 

Since D~ is a countable union of sets of finite one-dimensional Lebesgue measure, it has 
two-dimensional Lebesgue measure zero. However, its closure D~ may have a positive 
measure. Indeed that is the case for our parabolic map, and the following proposition is 
true. 

Proposition 1. For a parabolic map M with discontinuity, the set D~ is dense in the 
torus T 2 and its closure D~ has the full measure. 

Proof. Let B Xy£ be the e-ball around any point x on the torus. We have to show that 
it must intersect with D~ . Suppose that this e-ball and all of its images M n {B X)E ) do not 
intersect with D. Due to the shearing effect of the map the diameter of M n {B x ^ £ ) grows 
with n as 

diam(M n ( J B Xi£ )) = 2eJn 2 A 2 (a + ^ + 1. (12) 

Therefore if n is large enough, the diameter will be larger than \[2 and the image has to 
intersect the boundary D. 

Figure |5| shows the set D~ for the irrational map M^i defined by a = A = computed 
by a method similar to that used in |[(]], (the map M^i is distinguished by the fact that 
one element of the matrix g is equal to 7). Proposition 1 allows us to formulate main 
result of the present paper concerning dynamics of a generic parabolic map. 

Corollary. The dynamics generated by a generic parabolic map M has sensitive 
dependence on initial conditions in the torus T 2 . 
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This property is understood as in the study of Ashwin M regarding elliptic discontin- 
uous maps on the torus. 

Sensitive dependence on initial conditions in a set means that for all x\ £ Q there 
exists an e > such that for all 5 > there exists a point X2 € Q with \x± — a?2 1 < $ and 
|M n (xi) — M n {x2)\ > e for some n. 

Proof. As in the proof of the Proposition 1, any e-ball will intersect D after some 
iterations. Then the next iteration will split the ball into two disjoint sets separated by A, 
which is positive for a generic map M. Hence it has sensitivity on the initial conditions. 

Note that the proofs do not require that the map be irrational. It seems, however, 
that in the case of irrational maps, there is also another mechanism for the sensitivity 
due to the motion on the lines parallel to the eigenvector, whereas the rational maps do 
not have sensitivity along the eigenvector. Moreover, the above reasoning shows that the 
discontinuous parabolic maps fulfil a slightly different property of sensitive dependence 



introduced in the paper of Guckenheimer [12|. In any case, we emphasize that for this 



parabolic systems the sensitivity on initial conditions is not related with the Lyapunov 
exponents (which are zero in this case), but due to the discontinuity of the map. 

6 Diffusion and the decay of correlations 

Decay of correlations for a family of linear invertible parabolic maps on the torus was 
recently analyzed by Courbage and Hamdan |13|| . Although these systems are not chaotic 
(zero Kolmogorov-Sinai entropy), the correlations decay fast. While systems with sub- 
exponential correlations decay are found to be generic, a class of systems with exponential 
decay rate was found. 

Analogously, the non-invertible parabolic maps, discussed throughout this paper, dis- 
play fast decay of correlations, for generic, irrational maps. Figure ^(a) presents the 
normalized q~q autocorrelation function 

q{} ~ ((* - (%m ' (13) 

where the average is taken over initial points with respect to the natural measure. Numer- 
ical data are obtained from an ensemble of 10 6 initial points chosen randomly according 
to the natural measure /z* for the irrational map M 7 2. The fast decay of correlations, 
visible in this plot, was also observed for several other initial conditions. This suggests 
that the autocorrelation function is independent of the initial conditions, which supports 
the conjecture that the dynamics is ergodic with respect to the invariant measure As 
before, the fast correlation decay cannot be attributed to the positive Lyapunov exponent, 
but rather to discontinuity of the map M 7 2 . 

Figure ||(b) shows the Fourier transform of the autocorrelation function in (a), which 
is equal to the power spectrum of the g-time series. Notice the broad band continuous 
component in the spectrum similar to what is normally seen for chaotic maps. Also notice 
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that there are several peaks, suggesting some periodic components both in g-time series 
and the q—q autocorrelation function. 

In order to discuss the diffusion process it is convenient to work with the unfolded 
dynamics on the plane (i.e. the operation modulo 1 in Eq. (||) is not performed and the 
variable P in eqn.(|6|) takes arbitrary real values). Diffusion coefficient may be defined as 

D„ = lim D„(n) 

_ lim «* -<*»■>. (14) 

n—>oo ft 

Since P n +i — P n = (A/a)q n the position-position correlation function C g is equal to A/ a 
times the force-force correlation function in the direction p. This function determines the 
diffusion rate D p along the p axis. An argument similar to one in [|ll]] shows 

1 n 

D p = lim - £ {(A Pi - (A Pl )){A Pj - (A Pj ))) 

n >oo fi 

1,3=1 

{Api := P - Pi_i) 

-y n n—i 



}™L -E E ((A Pi -<Ap i ))(Ap i+i -(Ap j+i ») 



i=l j=l—i 
1 n n ~i /4\2 
i^-E E - ((QO - (qo)) 2 )C q (j) 
i=l 3=1— i 

-) 2 ^ E f 1 - -) <(«) - (®» 2 )^(j) 



-) 

a J 



j=l-n 

2 oo 

) ((Qo -(9o)) 2 ) E C q (n), (15) 



if the above sum converges, i.e. if the correlations decay not slower than n~ 2 . 

To study the the diffusion rate, we took 10 5 initial points distributed according to 
the natural measure and iterated them n times by a parabolic map M, and computed 
D p (n) := ((P-(P)) 2 )/n. Figure | is the comparison of this and the sum of autocorrelation 
in (p!q ) performed for an irrational map. Figure ^ shows the variance a 2 := ((P — (P)) 2 } 
as a function of time for four different parabolic maps. In each case, the mean was taken 
over 10 5 initial points. Observe two different kinds of behavior: for irrational maps, the 
variance is proportional to time, a 2 ~ D p n, which corresponds to the normal diffusion. On 
the other hand, for rational maps we receive a 2 ~ -C^ra 2 , which is called ballistic diffusion. 
In this case, which is sometimes referred to as the anomalous diffusion, the relation (15) 



looses its sense: the right hand side of this equation does not converge and the limit in 
the definition of diffusion coefficient (14) does not exist. Examples of anomalous diffusion 



were reported long time ago for a particular choice of the parameters of the standard map 
and more recently, for the kicked Harper model 15]. 
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7 Concluding remarks 



Although properties of linear area-preserving parabolic maps on the plane seem to be 
well understood, considering the same maps on the torus (by imposing periodic boundary 
conditions) makes the behavior of the system more complicated. Generically, such systems 
are discontinuous and non-invertible, and those properties lead to relevant dynamical 
implications. 

For rational maps, a single trajectory is contained in a finite number of lines on the 
torus. This leads to the anomalous diffusion, if the associated map is considered on the 
plane. On the other hand, for a generic, irrational map, a typical single trajectory explores 
the entire maximal invariant set. The dynamics shows sensitivity on the initial conditions 
and the diffusion process is normal. 

It is worth noting that several questions concerning the parabolic maps remain open. 
In particular, it would be interesting to prove the conjecture posed in this work: the volume 
of the maximal invariant set is a nowhere continuous function of the system parameter 
a; it is positive for all rational values of this parameter and zero for any irrational a. 
Moreover, our hypothesis that the dynamics of an irrational map is ergodic with respect 
to the natural invariant measure /i* is based on numerical experiments and still awaits a 
rigorous proof. 
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(a) 

gure 2: (a) Invariant set S of the symmetric parabolic map M s ;(b) magnification of S 
the vicinity of the critical point Q. 
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Figure 3: Periodic ID map derived from the symmetric map M s for q = 5/7. 
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(a) 

Figure 4: The maximal invariant set S 1 for the irrational parabolic map M 1 defined by 
a = 7 (a); and the rational maps related to the Fibonacci sequence a = 5/8 (b); and 
a = 2/3 (c). 
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Figure 5: The set D~ of the pre-images of the discontinuity lines up to k = 50 iterates 
for the system M^i. Thick lines correspond to the pre-images M~ 2 k (D) for low values 
of k, but the thickness is merely an artifact of the numerical procedure. The picture is 
produced using the 800 x 800 mesh. 
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Figure 6: (a) Force- force correlation function C q (n) for the system M 7 2 (b) and its Fourier 
transform S(lo). 
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Figure 8: Diffusion in parabolic maps on the torus, (a) golden map M^i (b) approximation 
to M y 2 with 7 = 3/5 (c) symmetric map M s (d) approximation to M s with a 
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(e) slope 1 (f)slope 2. The variance a grows linearly in time for the irrational map M^2 
and an irrational approximation of the symmetric map M s , while for rational maps its 
dependence is quadratic. 
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